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ABSTRACT 

The Constraint Network Analysis (CNA) web server 
provides a user-friendly interface to the CNA 
approach developed in our laboratory for linking 
results from rigidity analyses to biologically 
relevant characteristics of a biomolecular structure. 
The CNA web server provides a refined modeling of 
thermal unfolding simulations that considers the 
temperature dependence of hydrophobic tethers 
and computes a set of global and local indices for 
quantifying biomacromolecular stability. From the 
global indices, phase transition points are identified 
where the structure switches from a rigid to a floppy 
state; these phase transition points can be related to 
a protein's (thermo-)stability. Structural weak spots 
(unfolding nuclei) are automatically identified, too; 
this knowledge can be exploited in data-driven 
protein engineering. The local indices are useful 
in linking flexibility and function and to understand 
the impact of ligand binding on protein flexibility. The 
CNA web server robustly handles small-molecule 
ligands in general. To overcome issues of sensitivity 
with respect to the input structure, the CNA web 
server allows performing two ensemble-based 
variants of thermal unfolding simulations. The web 
server output is provided as raw data, plots and/or 
Jmol representations. The CNA web server, access- 
ible at http://cpclab.uni-duesseldorf.de/cna or 
http://www.cnanalysis.de, is free and open to all 
users with no login requirement. 

INTRODUCTION 

Proteins carry out their biological functions by interacting 
with other biomacromolecules or smaU molecules (1). These 



interactions require a certain degree of conformational 
adaptation to better complement the binding partners. 
That way, structural flexibility of proteins is linked to mo- 
lecular recognition and catalysis (2). Additionally, flexibil- 
ity (and its opposite, rigidity) also plays a central role for a 
protein's structural stabihty (3). Particularly, thermophilic 
proteins are in general more rigid than their mesophilic 
homologues to retain their fold at higher temperatures 
(4). Hence, knowing what can move in a protein is import- 
ant for hnking a protein's structure to its function and 
(thermo-)stabihty. Finally, information on protein flexibil- 
ity is increasingly incorporated in computer-aided drug dis- 
covery and design projects (5). 

X-ray crystallography, cryo-electron microscopy, single 
molecule fluorescence and nuclear magnetic resonance 
spectroscopy are experimental means from which the flexi- 
bility of a protein can be inferred (6-9). As an alternative, 
computational methods such as molecular dynamics 
(MD) simulations (10) and normal mode analysis are 
used to probe protein flexibility and dynamics (11). As yet 
another computational approach, a computationally 
highly efficient graph theory-based rigidity analysis for 
probing protein flexibility has been implemented in the 
Floppy Inclusions and Rigid Substructure Topography 
(FIRST) software (12). FIRST builds a constraint 
network from a biomolecular structure and then decom- 
poses this network into rigid clusters and flexible regions 
by using the pebble game algorithm (13,14). In the 
constraint network, atoms are represented as bodies, and 
covalent bonds and non-covalent interactions (including 
hydrogen bonds, salt bridges, stacked rings, and hydropho- 
bic tethers) are represented as bars (constraints) between 
them. Building on the ideas of Rader et cd. (15) on 
diluting non-covalent constraints in a constraint network, 
our group has developed the Constraint Network Analysis 
(CNA) approach that performs thermal unfolding simula- 
tions of proteins (16). CNA goes beyond the mere identifi- 
cation of rigid clusters and flexible regions in a biomolecule 
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in that it (i) provides a refined modeling of thermal unfold- 
ing simulations that also considers the temperature depend- 
ence of hydrophobic tethers; (ii) allows performing rigidity 
analyses on ensembles of network topologies, either 
generated from structural ensembles or by using the 
concept of fuzzy non-covalent constraints; and (iii) 
computes a set of global and local indices for quantifying 
biomacromolecular stability (17). Furthermore, CNA 
robustly handles small-molecule ligands in general. CNA 
has been successfully used for investigating protein 
thermostabihty, identifying unfolding nuclei ('structural 
weak spots'), and hnking protein flexibihty and function 
(18-20). 

In this study, we present the CNA web server that allows 
(i) setting up a variety of constraint network representa- 
tions of proteins from either single structures or ensembles 
of structures; (ii) performing rigidity analyses and thermal 
unfolding simulations on these networks; and (iii) process- 
ing the results. For that, the CNA web server provides a 
layer of user-friendly input and output interfaces around 
the CNA software. As input, the web server only requires a 
Protein Data Bank (PDB) code or user-provided PDB 
file(s) of the input structure(s), and choosing the simulation 
type. Results are presented in the browser in an interactive 
manner. For global indices, plots are provided; for local 
indices, plots and mappings onto the 3D structure [via a 
JmolApplet (http://jmol.sourceforge.net)] are provided. 
Weak spots predicted for a protein structure are also 
mapped onto the 3D structure. To the best of our know- 
ledge, there are no other web servers that allow performing 
and analyzing thermal unfolding simulations of proteins in 
as much detail as the CNA web server does. Of the two most 
closely related web servers, KlNARl (http://kinari.cs. 
umass.edu/Site/kinariWeb.html) (21) only performs rigid 
cluster decompositions, but no thermal unfolding simula- 
tions. Flexweb (http://flexweb.asu.edu/software/first/) does 
allow performing thermal unfolding simulations; however, 
it neither supports the use of ensembles of structures (which 
leads to more robust results from rigidity analyses) nor does 
it automatically determine phase transition points ('melting 
points') and unfolding nuclei (which extends the applica- 
tion domain of rigidity analysis to data-driven protein 
engineering). Thus, in addition to characterizing the distri- 
bution of flexible and rigid regions in a protein, the CNA 
web server can be used for probing changes in the flexibihty/ 
(thermo-)stability of a protein due to mutations or on 
ligand binding, and it aids in identifying structural weak 
spots in a protein that, when mutated, may improve the 
protein's (thermo-)stabihty. Note that the CNA web 
server solely characterizes what can move in a protein but 
does not simulate actual protein movements. To do the 
latter, the user is referred to the NMSim web server 
(http://cpclab.uni-duesseldorf.de/nmsim/ or http://www. 
nmsim.de) developed by us (22). 

MATERIALS AND METHODS 

Constructing a constraint network 

Proteins are modeled as body-and-bar networks where 
each atom is represented as a body with six degrees of 



freedom (23). Interactions between the atoms (covalent 
and non-covalent bonds) are modeled as a set of bars 
that restrict internal motion between the atoms. A 
covalent single bond is modeled with five bars allowing 
for the dihedral rotation about it; peptide and double 
bonds are modeled with six bars, disallowing any bond 
rotation. For example, a diatomic molecule with a single 
bond, owing to the five constraints, has seven degrees of 
freedom (6x2— 5), six of which represent the trivial 
overall rotations and translations and one of which repre- 
sents the internal rotation around the single bond. Non- 
covalent interactions, which contribute significantly to 
protein stability, are also modeled as bars. As such, 
hydrogen bonds (and salt bridges) are modeled with five 
bars, whereas hydrophobic and ring stacking interactions 
are modeled with two and three bars, respectively (24). 
Figure lA shows the structure of thermoly sin-like 
protease (TLP; PDB code: INPC) from which a body- 
and-bar network is then generated (Figure IB). 

Performing a rigid cluster decomposition 

Once the constraint network of a protein is built, the pebble 
game algorithm (13,14) as implemented in the FIRST 
software (12) decomposes the network into rigid clusters 
and flexible regions (Figure IC). The pebble game algo- 
rithm then computes the rigidity of the protein network 
at a bond level by determining whether a bond is part of 
a rigid cluster or a flexible region. As such, a rigid cluster is a 
set of atoms for which no internal motions are allowed and 
that move together in a collective manner. 

Simulating thermal unfolding 

By successively removing non-covalent constraints from a 
biomolecular constraint network, new network represen- 
tations at elevated temperatures are constructed. To do so, 
for a given network state s = f(7), hydrogen bonds 
(including salt bridges) with an energy ii/,/, > E^.,,, j,/, are 
removed from the network (25). This follows the idea 
that stronger hydrogen bonds wiU break at higher tem- 
peratures than weaker ones. The hydrogen bond energy 
scale Eh/, (25) is converted into a temperature scale T, 
using a linear equation proposed by Radestock and 
Gohlke (19). By default, the number of hydrophobic 
tethers is kept constant during the thermal unfolding 
simulation. However, as the strength of hydrophobic 
interactions increases with increasing temperature 
(26,27), hydrophobic tethers can also be treated in a tem- 
perature-dependent manner on request (20). Finally, a 
rigid cluster decomposition is performed on each con- 
straint network state .y to compute global and local flexi- 
bihty indices (Figure IC). 

Computing global and local flexibility indices 

The CNA web server computes in total six global and three 
local flexibility indices. Detailed definitions of the global 
and local flexibility indices are given elsewhere (17). In 
short, the global indices represent the macroscopic 
network flexibihty (and rigidity) of each network state 
,v: (i) The floppy mode density <P refers to the number of 
internal independent degrees of freedom that are associated 
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Figure 1. Covalent and non-covalent interactions in a protein structure (A) are modeled as bars in a body-and-bar network (B). A rigid cluster 
decomposition is carried out for all network states during a thermal unfolding simulation (C) and then post-processed to calculate flexibility indices, 
phase transitions, and weak spots. (D) Submission page to the CNA web server. 



with dihedral rotations, normahzed by the number of 
overall internal degrees of freedom associated with the 
total number of bodies in the network; (ii) the mean rigid 
cluster size S is computed with the size of the largest rigid 
cluster always being excluded; (iii) & (iv) the rigidity order 
parameter denotes the fraction of the network belong- 
ing to the giant percolating cluster (type 1) or the actual 
largest rigid cluster (type 2). The giant percolating cluster is 
the largest rigid cluster in the network state at the lowest 
temperature, i.e. with all constraints in place. During the 
thermal unfolding simulation, the melting of the giant 
percolating cluster is monitored, and the largest rigid 
subcluster of the previous giant percolating cluster 
becomes the new giant percolating cluster of the present 
network state s. In contrast, the actual largest rigid 



cluster is the largest cluster present at a network state s, 
irrespective of its evolutionary history during the thermal 
unfolding simulation; (v) & (vi) the cluster configuration 
entropy H (type 1 and type 2) is a measure of the degree of 
disorder in the realization of a given network state. 

Local indices characterize the flexibility of network at 
the bond level by monitoring the change in the flexibility 
for each bond during the thermal unfolding simulation: (i) 
the percolation index /?,■ is determined for each covalent 
bond by the E^^tMh value during the thermal unfolding 
simulation at which the bond segregates from the giant 
percolating cluster; (ii) the rigidity index is determined 
for each covalent bond by the E^utMh value during the 
thermal unfolding simulation at which the bond changes 
from rigid to flexible. For a atom-based representation 
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of a protein structure, the lower of the />, values (the 
average of the r,- values) of the two backbone bonds is 
taken; (iii) a stability map rcy is a 2D itemization of the 
rigidity index r,- and is derived by identifying the Ecui.hh 
value during the thermal unfolding simulation at which a 
rigid contact between a pair of residues represented by 
their Ca atoms is lost. Two residues are in rigid contact 
if they are part of the same rigid cluster. 

Identifying phase transitions and unfolding nuclei 

The global flexibihty indices are used for identifying phase 
transitions during the thermal unfolding simulation 
when the network switches from being largely rigid to 
flexible. Such transitions ('melting points') can be related 
to the (thermo-)stabihty of proteins (18-20). The CNA web 
server provides phase transition points for four of the global 
indices: Poo, type i and 2 and //type 1 and 2- A smoothcd spHne 
fitted to the global indices is used to identify the phase tran- 
sition points except for //type2. for which a double sigmoid 
curve is fitted. The phase transition points are identified 
as the maximum in the first derivative in the case of the 
smoothed spline and the maximum of the differences in 
the asymptote pairs in the case of the double sigmoid 
curve. The phase transition points can be further exploited 
to identify structural weak spots in the network from where 
an unfolding begins: these are residues that are part of the 
largest rigid cluster before the phase transition and become 
flexible afterwards. Weak spots provide hints as to where 
introducing mutations in a structure may improve a 
protein's (thermo-)stabihty. 

Simulating thermal unfolding of a network ensemble 

Results from rigidity analyses are in general sensitive to the 
input structure in that small changes in the conformation 
can lead to a different rigid cluster decomposition (28,29). 
To overcome this drawback, an ensemble-based variant of 
CNA has been developed that makes use of an ensemble of 
structures derived from MD simulations (20,28). Structural 
ensembles from any other sources can be used as input, too. 
As yet another alternative, the CNA web server provides 
an option to create an ensemble of networks from a sin- 
gle input structure by using fuzzy non-covalent constraints 
(Pfleger, C, Gohlke, H., unpublished data). Here, the 
number and distribution of non-covalent constraints are 
modulated by random components within certain ranges, 
that way simulating thermal fluctuations of a biomolecule 
without actually moving atoms. This approach avoids 
the use of computationally expensive MD simulations. 
For ensemble-based CNA, averages are computed over 
the entire ensemble for phase transition points and local 
indices. In the case of weak spots, the frequencies of occur- 
rence across the entire ensemble are reported. 

DESCRIPTION OF THE WEB SERVER 

Input 

The submission page to the CNA web server is shown in 
Figure ID. The CNA web server requires either a single 
structure of a protein provided as a PDB file or a PDB-ID 
(in which case the PDB file will be downloaded from the 



Research Collaboratory for Structural Bioinformatics re- 
pository), or multiple protein structures provided as PDB 
files in a compressed folder (allowed file formats are *.tgz, 
*.tar.gz or *.zip) or as a 'multi-PDB file' using MODEL/ 
ENDMDL cards. Furthermore, the web server requires 
selection of an analysis type and input of a given 
security code to prevent misuse. Analyses can be per- 
formed on a single constraint network derived from a 
single structure, on an ensemble of networks derived 
from a single structure or on an ensemble of networks 
derived from a structural ensemble. According to the 
selected analysis type, default parameters will be 
provided that have been successfully used in previous 
studies (17-20). In addition, the user can request the 
time-consuming computation of a stability map. An 
email address can be provided in which case a hnk to 
the results will be sent to that address. In either case, a 
hnk to a results page is provided after job submission for 
monitoring the progress of the computations and viewing 
the results in the web browser. The results will also be 
stored on the server for 10 days and can be accessed via 
the provided Hnk. 

If the user provides a PDB ID, missing sidechain and 
hydrogen atoms will be automatically added by the leap 
program of the AMBER suite (30). For amino acid 
sidechains, a standard protonation state is assumed, i.e. 
Asp and Glu are treated as deprotonated, and Arg and 
Lys as protonated. By default. His is singly protonated 
with the hydrogen on the epsilon nitrogen. Alternate 
sidechain conformations are allowed, but only the first 
conformation is considered. Water molecules and ligands 
will be removed. 

To have complete control on the input structure(s), a 
user can upload a/multiple PDB file(s). In that case, water 
molecules and ligands are considered, and correct bond 
orders for hgands are identified automatically using the 
Open Babel wrapper (Pybel) for Python (31,32). The 
correct identification of hgand bond orders is important 
to accurately set up covalent constraints for the network. 
In this case (i) hydrogen atoms must have been added 
already, (ii) alternate side-chain conformations are not 
allowed; and (hi) the protein structure must be complete, 
i.e. no missing sidechains are allowed. Finally, if the 
analysis is based on an ensemble of structures, aU 
protein structures must be identical except for the coord- 
inates. Otherwise, an error message will be issued. 

Output and representation of results 

A typical analysis run of the CNA web server for a con- 
straint network of a single structure takes a few minutes; 
for the analysis of an ensemble, several hours of 
computing time may be required. After the start, the pro- 
gression of the analysis is reported on the results page. On 
completion of a job, the results are presented in the web 
browser and, if an address is provided by the user, an 
email is sent with a link to the results, too. 

The first part of the results page contains a summary of 
input parameters that were provided by the user and a 
download link to the modified input file. The modified 
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input file contains all changes made during preparation of 
the uploaded structure. 

The second part of the results page contains a table that 
provides information about identified phase transition 
points. In the case of a single-network analysis, six plots 
of the global flexibihty indices are depicted: floppy 
mode density fp (Figure 2A), mean rigid cluster size 
5' (Figure 2B), rigidity order parameter P^o. type 1 and 2 
(Figure 2C and D) and cluster configuration entropy 
-fftype 1 and 2 (Figurc 2E and F). Identified phase transition 
points are marked by red vertical lines. In the case of an 
ensemble-based analysis, a summary of the identified 
phase transition points is presented instead. In any 



case, a file with the raw data can be downloaded by a 
Hnk given next to the headhne of the chapter for further 
evaluation. 

The third part contains two plots of the local flexibility 
indices for each residue: percolation index (Figure 3A) 
and rigidity index r, (Figure 3B). In the case of an 
ensemble-based analysis, the standard errors are 
depicted in gray in addition to the mean values. 
Furthermore, both indices are mapped onto the input 
structure in a color-coded fashion and shown in 
JmolApplets. A file with the raw data can be downloaded 
by a hnk given next to the headline of the chapter for 
further evaluation as can be two PDB files. These PDB 




Q. 

2^ 



-4 -3 -2 
E(,ut [kcal mol'^] 




(U 
Q. 



-4 -3 -2 
Ecut [kcal mol "^] 




[kcal mor^] 



[kcal mol "^] 



03 




[kcal mor^] 



Q. 




[kcal mol'^] 



Figure 2. Global indices for the thermal unfolding of TLP as a function of the hydrogen bonding energy cutoff £„„,/,/,: (A) floppy mode density <f>; 
(B) mean rigid cluster size 5; rigidity order parameter (C) type 1 and (D) type 2; cluster configuration entropy H (E) type 1 and (F) type 2. The 
red vertical lines (C-F) correspond to the identified phase transitions. 
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Figure 3. (A) Percolation index p, for TLP. Tlie lower the longer is a residue part of the giant percolating cluster during the thermal unfolding 
simulation. (B) Rigidity index /-, for TLP. The lower the the longer is a residue part of a rigid cluster during the thermal unfolding simulation. 
Red- and green-dashed horizontal lines represent the identified phase transition point and the working temperature of TLP, respectively. The central 
ff-helix and two preceding Gly residues (residues 136-154) residues are enclosed in a red rectangle. On the right, the respective indices are mapped 
onto the input structure in a color-coded manner. (C) Stability map cc,, for TLP. Red colors indicate pairs of residues where no or only a weak rigid 
contact exists. In contrast, blue colors indicate strong rigid contacts. The black box with a continuous line covers the N-terminal giant rigid cluster, 
whereas the box with the broken line indicates a rigid cluster in the C-terminal domain. (D) Weak spots in the TLP structure are represented by red 
spheres. 
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files contain the modified input structure of the protein 
with the respective index values of atoms in the 
B-factor column. 

The fourth part contains a stabihty map rc,y (Figure 3C) 
(if requested) and information about weak spots identified 
in the protein (Figure 3D). For rcjj, the color code reveals 
for each residue pair if the rigid contact is weak (red) or 
strong (blue). Information about weak spots is mapped 
onto the input structure and displayed in a JmolApplet: 
in the case of a single-network analysis, identified weak 
spots are marked by red spheres; in the case of an 
ensemble-based analysis, the frequency of a residue for 
being identifed as a weak spot across the entire ensemble 
is depicted. This information is provided as a plot, too. 
Again, two files can be downloaded by links given next to 
the headline of the chapter for further evaluation: a file 
with the raw data and a PDB file containing the modified 
input structure of the protein with the respective frequency 
values in the B-factor column. 

Implementation 

The CNA web server has been implemented in Python, as 
have been the underlying CNA routines (16). All plots are 
generated with Gnuplot. Given the low computational 
demand of our approach, up to 10 submitted jobs can 
be run in parallel at present. 

APPLICATION TO TLP AS A TEST CASE 

In a previous study from our group, the CNA approach 
has been shown to link protein rigidity and 
thermostability by correctly predicting which one of two 
proteins is more thermostable for two thirds of a data set 
of 19 pairs of proteins from mesophilic and thermophific 
organisms (18). In the same study, we showed that struc- 
tural weak spots identified by CNA agreed with the pos- 
itions of thermostabilizing mutations for two systems 
investigated. Extending the application domain of CNA 
one step further, the local flexibihty and rigidity distribu- 
tions in the active sites of 3-isopropylmalate dehydrogen- 
ase and TLP were then linked to enzymatic activity at 
different temperatures (19). Recently, using CNA on 
citrate synthases with distinct thermostabilities, we 
showed that thermostabilizing mutations preferentially 
occur at structural weak spots with a high mutation 
ratio (14). 

Here, we demonstrate the appHcation of the CNA web 
server to TLP from Bacillus cereus (PDB code: INPC) 
using a single-network approach with default settings. 
These results are also available as a sample run on the 
web server. Plots for global indices are shown in Figure 
2. The decrease in network rigidity of TLP with increasing 
temperature (equivalent to a decreasing is evident 

from the plots of all six global indices. Interestingly, the 
decay of the giant rigid cluster occurs in a hierarchical 
fashion as reflected by the presence of multiple steps in 
the Poo and H profiles (Figure 2C-F). The reason for this 
percolation behavior is that the TLP structure is 
composed of multiple sub-domains (N-terniinal j6-sheet 
domain and C-terminal a-hehcal domain connected by a 



central helix) that segregate from the giant cluster inde- 
pendently from each other. Phase transition points at 
which the TLP structure sharply loses rigidity on 
thermal unfolding are then identified: an early phase tran- 
sition point at E,,i„i,i, = — 1.14kcalmol~' is identified from 
the /"oo and //typei profiles when the C-terminal a-helical 
domain segregates from the giant cluster; in contrast, 
Htyps 2 fosters the identification of a late transition at 
= — 2.55 kcalmoP' [or, equivalently, 351 K (19)], 
which represents the final substantial decay of the rigid 
core. Such late transitions have been found to best relate 
to the melting of a protein. The temperatures of the tran- 
sition points should be considered relative values only 
because the absolute values may depend on the size and 
architecture of the analyzed protein (19). Still, the tem- 
peratures are helpful, e.g. when it comes to comparing 
the thermostabihty of two or more homologous proteins 
or the stabihty of a wild-type with its mutant (18-20). 
Accordingly, the //type 2 phase transition point for a 
thermophilic homologue of TLP, thermolysin from 
Bacillus thermoproteolyticiis, was found at 373 K (24), 
~20K higher than for the mesophilic homologue 
investigated here. 

The percolation index plot depicts the percolation 
behavior of TLP at a residue level (Figure 3A). The 
C-terminal a-helical domain segregates from the giant 
cluster early during the thermal unfolding simulation at 
/i™/,/i/> ~ — 1-1 kcalmol"'. The giant percolating cluster, 
which then consists of mainly the /J-sheet region and an 
a-helix in the N-terminal domain, disintegrates into 
smaller clusters at the later phase transition point 
(/^(•i//./;/; = — 2.55 kcalmol"'). Unsurprisingly, the rigidity 
index, which monitors the flexibihty at a residue level, 
shows secondary structure elements to be more rigid than 
loops. As an exception, residues 1 17-120 forming a loop in 
the N-terminal domain are always part of a rigid cluster 
throughout the thermal unfolding simulation owing to a 
network of hydrophobic tethers. For TLPs, the central a- 
helix (residues 139-154) and the preceding Glyl36 and 
Glyl37 are important with respect to a postulated hinge 
bending motion (33,34), suggesting that these residues 
should be flexible at the working temperature of TLP 
[342 K, equivalent to = -2.1 kcalmol"' (19)]. 

This is also found in the rigidity index profile (Figure 
3B). Furthermore, contacts of these residues with other 
residues of TLP are less stable than contacts between 
residues of the giant cluster (black box with continuous 
fine) and contacts between residues of a large rigid cluster 
in the C-terminal domain (black box with broken fine) as 
identified by the stability map (Figure 3C). This is again in 
line with the hinge character of the central a-helix. Finally, 
several of the weak spots identified by CNA in the N- 
terminal j6-sheet domain of TLP (Figure 3D) have been 
shown to improve the protein's thermostabihty on 
mutation in previous studies (35-38). 

CONCLUSIONS 

Increasing evidence of the importance of protein flexibility 
has warranted the development of efficient and accurate 
computational tools for characterizing protein flexibility 
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at global and local levels. To this end, the CNA approach 
has been developed for deriving maximal advantage from 
information on biomolecular flexibility and rigidity by 
linking results from rigidity analyses to biologically 
relevant characteristics of a structure, such as (thermo-) 
stabihty and function (16). In particular, CNA provides a 
refined modehng of thermal unfolding simulations that 
also considers the temperature dependence of hydropho- 
bic tethers, allows performing rigidity analyses on ensem- 
bles of network topologies and computes a set of global 
and local indices for quantifying biomacromolecular sta- 
bility. Furthermore, CNA robustly handles small- 
molecule hgands in general. To make these computations 
available for users even with only minimal or no prior 
knowledge of structural bioinformatics techniques, we de- 
veloped the CNA web server. It provides a user-friendly 
interface, requires minimal input and displays the results 
intuitively both as plots and mappings onto the protein 
structure via JmolApplets. As a typical analysis run of the 
CNA web server for a constraint network of a single struc- 
ture only takes a few minutes, we strongly believe that the 
CNA web server will be a valuable (interactive) tool for 
data-driven protein engineering and estimating the influ- 
ence of hgand molecules on biomolecular stability. 
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